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ABSTRACT 

In two long-duration balloon flights in 2003 and 2006, the TRACER cosmic- 
ray detector has measured the energy spectra and the absolute intensities of the 
cosmic-ray nuclei from boron [Z = 5) to iron [Z = 26) up to very high energies. 
In particular, the second flight has led to results on the energy spectrum of the 
secondary boron nuclei, and on the boron abundance relative to that of the heav- 
ier primary parent nuclei, commonly quantified as the "B/C abundance ratio". 
The energy dependence of this ratio, now available up to about 2 TeV per amu, 
provides a measure for the energy dependence of cosmic-ray propagation through 
the Galaxy, and for the shape of the cosmic-ray source energy spectrum. We use 
a Leaky-Box approximation of cosmic-ray propagation to obtain constraints on 
the relevant parameters on the basis of the results of TRACER and of other mea- 
surements. This analysis suggests that the source energy spectrum is a relatively 
soft power law in energy E~ a , with spectral exponent a = 2.37 ± 0.12, and that 
the propagation path length A(E) is described by a power law in energy with 
exponent 5 = 0.53 ± 0.06, but may assume a constant residual value A at high 
energy. The value of Ao is not well constrained but should be less than about 
0.8 g cm -2 . Finally, we compare the data with numerical solutions of a diffusive 
reacceleration model, which also indicates a soft source spectrum. 
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Introduction 



The composition and energy spectra of Galactic cosmic rays change while the particles 
propagate from the acceleration site to the observer. Thus, the cosmic-ray population ob- 
served in measurements is different from the source population, even in the approximation 
that the ambient population is the same in time and space throughout the volume of the 
Galaxy. To determine the characteristics of the cosmic-ray sources, the propagation effects 
must be deconvoluted from the measured data. 

Processes affecting the cosmic rays after they are released from their sources include 
diffusion and diffusive or convective escape from the Galaxy; interactions with the compo- 
nents of the interstellar gas leading to the loss of the primary and to the production of 
secondary particles; decay of radio-active components; secondary acceleration in the inter- 
stellar medium (ISM); ionization energy loss; or radiative energy loss (for electrons). Before 
the cosmic rays reach an observer near the Earth, they may also be affected by solar mod- 
ulation. Many of these processes (including their energy dependence) are more difficult to 
quantify at relatively low energies, below a few GeV per amu, than in the highly relativistic 
region. Nevertheless, even at high energies, the numbe r of parameters t hat appear in the 
full transport equation for Galactic cosmic rays (see e.g. iGinzburgl (119641 )) is too large, and 
the parameters are often too poorly known, as to permit an easy analytic description of 
the effects of cosmic-ray propagation. In the following, we will review recent experimental 
data in the context of the long-popular "Leaky-Box" approximation of particle transport in 
the Galaxy. We also will compare t he data with a diffusive propagation model using th e 
prescriptions of the GALPROP code ( jStrong and Moskalenkolll998l ; IVladimirov et al.ll2010l ). 



At low energies, observational data of good precision have been available for some time 
from detailed measurements with a number of spacecraft. For the high energy region, from 
GeV amu" 1 to TeV amu" 1 energie^], few measurements in space are available, but a number 
of long-duration balloon (LDB) flights in recent years have provided new data. The following 
study will therefore concentrate on the high-energy region. 

Specifically, this study is motivated by recent results from the TRACER cosmic-ray 
detector. This instrument is currently the lar gest and most sensitive detector for the heavier 
(Z > 3) cosmic-ray nuclei ( lAve et al.l 1201 ll ). TRACER was flown in two LDB flights, in 
2003 in Antarctica (LDB1), and in 2006 from Kiruna, Sweden (LDB2). The 2003 flight 



1 The measured quantity for the TRACER instrument (and many others) is the Lorentz-factor 7 = E/mc 2 , 
which is equivalent to (kinetic energy per amu + 1) if the mass is expressed in amu. Within the uncertainties 
of most measurements, the kinetic energy per amu is identical to the more commonly used quantity kinetic 
energy per nucleon. 
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yielded a comprehensive set of the spectra of the major primary nuclei from oxygen {Z — 8) 
to iro n (Z = 26), cov ering the energy range from a few GeV per amu to several TeV per 
amu (lAve et al.ll2008l ). For the 2006 flight of TRACER, several significant upgrades of the 
instrument permitted the inclusion of the lig ht elements of boron (Z = 5), carbon [Z = 6) 
and nitrogen (Z = 7) in the measurement (jObermeier et al.l 1201 ll ). Particularly impor- 
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20081 . 120091 ). the current data set on the composition of high-energy cosmic-ray nuclei should 



lead to more stringent constraints on the propagation of high-energy cosmic rays than pre- 
viously possible. 



2. Parameters of Cosmic Ray Propagation 

The "Leaky-Box" approximation parametrizes the diffusion of cosmic rays in the Galaxy 
by introducing the average propagation path length A(E), which is inversely proportional to 
the diffusion coefficient and may depend on energy E: 

A(E)=f3cpr(E), (1) 

where t(E) is the average containment time in the Galaxy, = v/c is the particle's velocity, 
and p is the mass density of the interstellar gas. Similarly, the average spallation path length 
A s (v4) quantifies the interaction of a cosmic-ray nucleus with mass number A by spallation 
reactions with interstellar gas nuclei of mass M: 

A.(A) = ^. (2) 

For relativistic energies, we assume that the spallation cross section a depends only on A, 
ignoring a possible small energy dependence. The values of A and A s are of the same order 
of magnitude at GeV energies. 

For highly relativistic nuclei, one may disregard ionization energy losses and solar mod- 
ulation. If one also ignores convective motion of particles as well as secondary acceleration 
and radio-active decay, one obtains a conti nuity equat i on for the ambient differential density 



Ni of cosmic-ray species i in the form (e.g. lAve et al.l (120091 )): 
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Here, Qi(E) is the rate of production in the source, and \-i is the differential spallation 
path length for a nucleus k to spallate into i. As is common practice, we will accept this 
form of the Leaky-Box equation as a valid approximation for the energy spectra of the more 
abundant cosmic-ray species. 

The Leaky-Box approximation assumes a continuous distribution of cosmic-ray sources 
and makes no assumption about the boundaries of the d iffusion region and t he extent of 



the Galactic halo. Nevertheless, it has been found (e.g. iPtuskin et al.l (120091 )) that more 
realistic diffusion mod els, such as the numerical in tegration of the transport equation in 
the GALPROP code (jStrong and Moskalenkd Il998l ). lead to results for the major stable 



cosmic-ray nuclei, which are equivalent to the Leaky-Box predictions at high energy. 

The propagation path length A decreases with energy above a few GeV amu" 1 . The 
first evidence for this phenomenon came from the observation of the dec reasing relative 



abundance of secondary, spallation-produced cosmic-ray nuclei with energy (IJuliusson et al. 



1972 



Smith et al.l I1973T). This and s ubsequent measure ments, in particular on HEAO 



3 (jEngelmann et al.l 119901 ) and CRN (jSwordy et al.l Il990l ) , indicated that the energy de- 
pendence of A might have the form of a power law E~ 5 , with 8 ~ 0.6, at relativistic energies 
(> 20 GeV amu' 1 ). 

It is of historic interest that at about the same time, stochastic shock-acceleration in 
supe rnova- remn ants was proposed as an efficient process to produce Galactic cosmic rays 
(e.g. iBelll (119781 )). However, for strong shocks, this theory predicts in, to first order, a source 
energy spectrum Qi(E) in form of a power law E~ a , with a ~ 2, much harder than the 
well-known E" 2 - 7 behavior of the observed spectrum. A cosmic-ray escape from the Galaxy 
that scales with E~~ 6 indeed seemed to provide an easy solution to this dilemma. This fact 
may well have fostered the acceptance of the shock acceleration model. However, a closer 
look at currently available data as will be attempted in this paper, may reveal complications 
to this simple picture. 

The E~ s energy dependence of the propagation path length would lead to very small 
values of A at high energies that might not be consistent with the reported isotropy of the 
cosmic-ray flux. Therefore, one may speculate that the path length approaches a residual 
value A at high energies: 

A = C ■ E~ 5 + A . (4) 

Physically, the residual path length would characterize a minimum column density of matter 
that a cosmic-ray particle must traverse, even at very high energy. This matter could be 
located near the cosmic-ray source, or it could signify the minimum distance to a source. 



It is commonly assumed, and consistent, to first order, with the shock acceleration hy- 
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pothesis, that the energy spectrum of all primary cosmic-ray nuclei, has the same dependence 
on energy (or rigidity), typically in form of a power law with common index a: 

Q l (E)=n l -E~ a , (5) 

Equation ([3]) then takes the form 

"'< B > = (C • E-> + A.)- + AM)- 1 (6) 



This equation describes the connection at high energies between the measured cosmic- 
ray spectrum Ni(E) for every stable nuclear species i, and the power-law spectrum and 
relative intensity n, of that species at the source. Besides a normalization factor C, this 
relation uses just three free parameters, the source index a, the propagation index 5, and 
the residual path length A . 

The equation illustrates in a simple way how the competition between diffusive and 
spallation loss, expressed by the term 1/(A(E)~ 1 + A s (A)~ r ), is reflected in the shape of 
the measured energy spectra: the smaller of the two A-parameters will most strongly affect 
the resulting spectrum. Current data, which will be further discussed below, and which are 
summarized in Figure HI show that A(E) is smaller than A s for much of the energy region 
of concern. However, to completely ignore the influence of spallation on the spectral shape 
we must require that A(E) is smaller than A s by at least a factor of ten. For the nuclei 
of concern here, this will be the case for energies in the 100 - 1000 GeV amu -1 region if 
the residual path length A is zero (as is illustrated in Figure 0J. At these energies, most 
measured data are affected by significant statistical errors and do not strongly constrain the 
functional shape of the observed spectrum. Returning to Equation (jBj), we predict that the 
observed spectrum will be softer than the source spectrum E~ a , but that the often used 
description of the observed spectrum as a power law E~ r with index T = (a + 5) can only 
become valid at the highest energies, and if the residual path length is zero. At lower energies 
(and that is where the most accurate data are available), a power- law fit with a smaller value, 
r < (a + 5) would approximate the observed spectrum. 



3. The Data 

The measurements available from the two long-duration balloon flights of TRACER 
cover the energy spectra of ten cosmic-ray elements from boron (Z = 5) to iron [Z = 26). 
The data of the first flight, for the primary nuclei from oxygen (Z = 8) to iron, have 
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be en reported by lAve et all (120081 ). and the results of the second LDB flight are described 
by lObermeier et al.l ( 1201 ll ). Where overlap exists, the two data sets agree well with each 
other. The resulting energy spectra are shown in Figure [TJ Figure [T] also indicates a simple 
power law fit with index 2.65 ± 0.05 which was found to describe all primary sp ectra quite 
well above 20 GeV amu" 1 , without any significant change with charge number Z (lAve et al. 



2009|). 



The data of the second flight include the energy spectra for the light elements boron 
(Z = 5) and carbon (Z = 6) up to about 2 T eV amu -1 . The resulting boron-to-carbon 
(B/C) abundance ratio (lObermeier et al.ll201lf) is shown in Figure [2] This figure also in - 



1990) 



20101) 



eludes previous data from HEAO-3 (Engcl mann et al.l Il990l ). CRN (ISwordv et al 
ATIC JPanov et al.ll2007bh . CREAM (Uhn et aliboosh . and AMS-01 (Uguilar et a" 
with their reported statistical uncertainties. At the highest energies, all results are based 
on few events, and not all measurements then define the statistical uncertainty in the same 
way. The results on the B/C ratio from balloon flights include a correction for atmospheric 
production of boron which may become sizable at high energy. The level of the correction 
for the TRACER data, which has been subtracted from the ratio (at an average residua l 
atmosphere of 5.2 g cm" 2 ), is indicated as a dashed line in Figure [2] (jMuller et al.l 1201 lh . 
Details about this correction will be published separately. 

Also shown in Figure [2], is a prediction for the B/C ratio corresponding to an energy de- 
pendence of the escape path length A, essentially proportional t o E~ os . Specifically, we have 



chosen for this prediction the parametrization of lYanasak et al.l (120011 ). which was developed 
to include low-energy data from HEAO-3 and ACE-CRIS (below ~ 10 GeV amu" 1 ): 



A(R) [g cm" 



C/3 



{PR) 5 + (0.714 ■ /3R)" 1 - 4 



+ A 



o- 



(7) 



where R = pc/Ze is the particle's rigidity and = v/cis the particle's velocity. lYanasak et al. 
(120011 ) use a propagation index 5 = 0.58, a residual path length A = 0, and a normalization 



C = 26.7. It should be noted that at high energies Equation (j7]) is equivalent to Equation 



4. Constraints on Propagation Parameters 



We now will attempt to derive constraints on the propagation parameters from the 
measurements just described, i.e. constraints on the cosmic-ray source spectral index a, and 
on the energy-dependence of the propagation path length char acterized by the parameters 
6 and A . The present work continues the fitting procedures by lAve et al.l (120091 ) applied to 
the data of the first LDB flight of TRACER in 2003. 
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Referring the reader to the paper by lAve et al.l (120091 ). we recall that this first flight led 
to the energy spectra of the heavier p rimary cosmic-ra y nuclei, but did not include results 
on secondary elements such as boron (lAve et al.ll2008l ). Hence, no new information on the 
propagation path length was obtained, and for the fitting procedure fixed parameters 6 = 0.6 
and C = 26.7 were assumed. However, Ao and a were treated as free parameters. The most 
striking feature of the measured data was the common power-law appearance of all measured 
energy spectra from 20 GeV amu" 1 to several TeV amu -1 (with the same index of 2.65). This 
feature could only be reconciled with the prediction of Equation ([6]), if the energy spectrum 
at the source was fairly soft, with a probable value of the source index a between 2.3 and 2.45. 
This value is considerably larger than the first-order expectation of a ~ 2.0 for acceleration 
in strong shocks. The data from this flight did not place strong constraints on the residual 
path length A , and could not exclude a non-zero value for this parameter. With these fitting 
results, the measured energy spectra of the individual elements were extrapolat ed back to 



the sources, and the relative elemental source abundances m had been obtained (Ave et al. 



2009|). 



With the measurement of the energy spectrum of the secondary nucleus boron, and 
of the secondary /primary intensity ratio, i.e. the B/C ratio, in the second balloon flight in 
2006, we now attempt to derive further detail. We use Equation (j3J), which for boron does 
not contain a source term Qj. Introducing an effective path length X^b (see Equation (jUJ)), 
the B/C ratio can then be expressed as: 



Nb 
N c 



X 



A- 



Ai 1 ' 



(8) 



Here, we further assume that boron is produced only by spallation of carbon and oxygen, 
i.e. the contributions from the spallation of nitrogen (amounting to just ~ 3% of the boron 
intensity) and from nuclei with Z > 8 are ignored. Finally, we assume that there are no 
significant contributions to the intensities of carbon and oxygen from spallation of heavier 
nuclei. These assumptions seem to be justified by the dominant intensities of carbon and 
oxygen among the primary nuclei. The effective production path length for boron X^b 
includes both carbon and oxygen as parent nuclei: 



A 



Ib = + No/Nc ■ Xo 1 . 



O^B- 



(9) 



The ratio Nq/Nq refers to the intensity ratio of the parent nuclei oxygen and carbon on 
top of the atmosphere. This ratio can be taken as independent of ener g y, and is close to 



unity (lObermeier et al-lbOlll Miiller et al.lll99ll : Engelmann et al.l 11990 : lAhn et al.l 120081 ). 



The spallation path le ngth A r in Equation ( | 8l) is derived from a g eometrical parametrization 
of the cross sections ( iBradt fc Petera Il950l : IWestfall et al.lll979l ). and the production path 
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l engths A in Equation (Q are derived from partial cross sections determined by I Webber et al. 
(Il990bl ). Specifically, we use A B = 9.3 g cm -2 , and X^ B = 26.8 g cm" 2 (assuming the 



interstellar medium as a mixture of 90% hydrogen and 10% helium by number). 

The fitting function is then given with Equation (|SD, with the escape path length A as 
expressed in Equation (J7J). Compared to using the high-energy form of Equation (j3J), this 
has the advantage that data below ~10 GeV amu" 1 can be included in the fit. The only 
unknown quantity in Equation (jSJ) is the energy dependence of the propagation path length 
A with the parameters S and A . 

We have fitted the data on the B/C ratio versus energy as measured by TRACER to 
a variety of values for S and A . A probability contour map of the fitting results is shown 
in Figure EJ The best fit for the propagation index is 5 = 0.53 ± 0. 06 g cm" 2 , and is quite 



close to the value of 0.6 which was used in the previous analysis of lAve et al.l (120091 ). The 
best value for the residual path length, A = 0.3li ; 3 i g cm" 2 , is less well defined, and still a 
solution with Ao = cannot be excluded within the present accuracy of the TRACER data 
alone. The corresponding escape path length A together with its uncertainties is shown in 
FigureHJas a function of energy. The figure indicates that a cosmic-ray nucleus most probably 
traverses a column density of 2.5 ±0.9 g cm" 2 of matter at an energy of 50 GeV amu" 1 before 
escaping the Galaxy. At 1000 GeV amu" 1 , the path length will be between 1.6 g cm" 2 and 
0.28 g cm" 2 , with a best-fit value of 0.76 g cm" 2 . For comparison, the figure also indicates the 
energy-independent spallation path lengths for the primary elements carbon and iron. The 
result of the fitting procedure is shown in Figure [5] as a solid line. The fit to the TRACER 
data alone overshoots the low energy data of other measurements by about 10%-20%. 

To refine the fit we may attempt to use the total data set currently available for all 
reported B/C ratio measurements at high energy (see Figure [2]) in the fitting routine. The 
result for the propagation parameters of this analysis essentially agrees with the analysis of 
the TRACER data alone, but leads to values which are more tightly constrained: we now 
obtain 5 = 0.64 ± 0.02, and Ao = 0.7 ± 0.2 g cm" 2 . If this is correct, it would be the first 
evidence for a non-zero residual path length. However, we feel that this conclusion must be 
taken with caution as the different measurements may be affected by different systematic and 
statistical uncertainties, including uncertainties in the cross sections and the contribution of 
nitrogen. The result of the fit to all data is also illustrated in Figure [5] (dashed line). 

Turning our attention now toward the shape of the energy spectra of primary nuclei, we 
show in Figure[6]the energy spectrum for oxygen, Nq{E) (multiplied with E 2S5 ), as measured 
in the two TRACER flights. The prediction of Equation ([6]) is applied with the propagation 
parameters from TRACER as just described {5 = 0.53±0.06, A = 0.31±g;|f g cm" 2 ). When 
the source spectral index a is varied, a best fit for the source spectrum is obtained with 
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a = 2.37 ± 0.12. The shape of this source spectrum is indicated in Figure [6j As mentioned, 
the production of oxygen by spallation is ignored in the fitting procedure, and the spalla- 
tion path length of oxyge n is derived from a geometric al parameterizatio n (IBradt fc Peters 



19501 ; IWestfall et al.l Il979f ) . As in our earlier analysis (I Ave et al.l 120091 ) , this fit again in- 



dicates a fairly soft energy spectrum at the source. One may note that the corresponding 
energy spectrum expected to be observed at the Earth is not a straight power law. The 
curvature reflects the effect of the non-zero residual path length. A second deviation from 
a straight power-law behavior of the observed spectrum is expected at lower energy due to 
the competition of spallation and escape during galactic propagation, and eventually due to 
solar modulation. Unfortunately, the uncertainties in the data are too large as to determine 
whether these subtle curvature effects can be real. 



5. Comparison with a diffusive propagation model 



A more realistic approach to describe the cosmic-ray propagation in the Galaxy involves 
numerical so lutions of the transport equat i ons. The most common ly used tool is the GAL- 
PROP code (jstrong and Moskalenko 1998 ; Vladimirov et al. 2010 ). It has been successfully 
employed to provide a self-consistent description of most observations of cosmic rays includ- 
ing protons and anti-protons, heavier nuclei, electrons (but apparently not positrons), over 
a wide range of energies, and of diffuse Galactic gamma rays. 

We apply GALPROP with input paramete rs (e.g. scale height, cross sections, source 
abun dances) as suggested in several publications (jStrong and Moskalenkoll200ll ; Ide Nolfo et al 
20061 ) . We also accept a diffusive propagation model with an energy-dependent diffusion 
coefficient commensurate with Kolmogorov turbulence in the Galactic magnetic field, i.e., 
D proportional to E s , with 5 ~ 1/3. The model also permits reacceleration of cosmic 
rays in interstellar space. This model is widely used to describe cosmic ra diation in the 



Galaxy and the values of its para meters have been applied for some time (IPtuskin et al. 



2006|; 


Strong and Moskalenko 


2001 


) and have recently been confirmed again by 


Trotta et al. 


(2011 


)• 







In this model, no asymptotic value for the diffusion coefficient is invoked, so there is no 
equivalent to the residual path length A . For this specific model, the preferred value for 
the power-law index of the cosmic-ray source spectra is a = 2.34, very similar to the value 
resulting from our Leaky-Box fit described above. 



We compare the prediction of this model with the measurements of the B/C ratio, as 
shown in Figure [TJ Indeed, the model describes the data well, and closely constrains the 
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value of the propagation parameter to 5 ~ 0.34. This value is smaller than the parameter 
5 = 0.53 ± 0.06 obtained for the Leaky-Box fit. 

The choice of parameters in the diffusion model is not unique. While it is not the 
purpose of this study to review the various parameter selections that have been discussed in 
the literature, we emphasize that reacceleration affects the cosmic-ray energy spectra in the 
diffusive model, and contributes to the energy dependence of the B/C ratio. If reacceleration 
were insignificant, the diffusion model would be commensurate with a value S ps 0.6 obtained 
in the Leaky Box ap proximation, but the source spectral index a would be smaller than 



before, a » 2.15 (e.g. IPtuskin et all (12006( 1 ). 



At the very highest energies, a non-zero value of the residual path length A would 
flatten the energy dependence of the B/C ratio. While this would be a distinct signature 
not predicted in the diffusion model (unless a primary contribution to boron is invoked), its 
observation is still within the fringes of current experimental accuracy. 



Discussion and Conclusion 



We have discussed constraints on the propagation of cosmic rays through the Galaxy 
derived from the measurements with the TRACER instrument, and from several other recent 
observations. A simple Leaky-Box model is used that reduces the number of free parameters 
to the spectral index a of the source energy spectrum (assumed to have the same power- 
law form for all primary nuclei), and to the energy dependence of the diffusive escape path 
length A esc , which decreases with energy as a power law with index S, but may also exhibit 
a constant residual value A (see Equation (j3J). 

A previous analysis had been performed based on the data on the primary cosmic- 
ray nuclei from oxygen (Z = 8) to iron (Z = 26), obtained with the first LDB flight of 
TRACER (lAve et al.l 120091 ). This analysis had led to the conclusion that the energy spectra 
of the cosmic rays at the sources must be quite soft, with a = 2.3 — 2.4. No new data on 
secondary cosmic rays were available at that time; hence, a propagation index S = 0.6 was 
assumed in the analysis, and no strong constraint on A could be obtained. 

The d ata discussed in the pr esent work include results from the second LDB flight of 
TRACER ( lObermeier et al.l 120111 ). which provides a new measurement of boron and carbon 
nuclei, and of the B/C ratio at high energies. Therefore, the fitting routine could include 
the parameters, a, 5, and Ao as free parameters. Remarkably, the best-fit values of these 
parameters agree well with the conclusions of the previous work: again, the source spectrum 
is quite soft, with a = 2.37 ± 0.12 and the propagation index 5 = 0.53 ± 0.06 is close to the 
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previously fixed value of 0.6. The residual path length Ao may be about 0.3 g cm 2 but still 
has a large uncertainty which cannot exclude a value A = 0. 

The agreement with the earlier work has an important consequence: the relative abun- 
dances of the individual elements at the sources (i.e., th e numbers rij in Equations ([5]) and 



([6])), which had been calculated in the earlier analysis of I Ave et al.l ( 120091 ). remain valid. In 



particular, the comparison o f these values with the sol a r system abund ances (often also called 



universal abundance scale ( jAnders fc Grevessd Il989l ; iLodderd 120031 )). and the correlations 



with parameters such as the first ionization potential, or the condensation temperature, do 
not need to be updated. 

The softness of the energy spectra at the source predicted by the present analysis agrees 
well with the source spectrum preferred by the diffusive reacceleration model. An open ques- 
tion is the significance of the residual path length. This question cannot be fully answered 
until a new generation of instruments provides measurements of the B/C ratio with greatly 
improved statistical accuracy in the TeV amu -1 region. 

The Leaky-Box model of the present analysis does not include reacceleration in inter- 
stellar space as a significant contributor among the effects of Galactic propagation. Con- 
sequently, the energy dependence of the propagation path length is fairly strong, with the 
index 5 close to the often used value of 0.6. 

Finally, one may consider deviations from a pure power law behavior of the cosmic-ray 
energy spectra either at the sources, or at the observation site. In the present analysis, a strict 
power law was accepted for the source spectrum. As Figure [6] indicates, it would then be 
inevitable that the observed spectrum deviates from a power law form, but the deviation is so 
small that it cannot be observed within current observational uncertainties. Or, if detected, 
it would be difficult to ascribe small effects in the observed spectra to either propagation 
or source effects. Once again, much improved statistical and systematic accuracy in future 
measurements is required to settle this issue. 

In summary, it appears that the recent cosmic-ray measurements at high energy are 
beginning to probe details of the source and propagation characteristics that have remained 
unexplored for a long time, and one may expect more answers to the remaining questions 
from a next generation of instruments. 
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Fig. 1. — Compilation of the differential energy s pectra mea s ured by T RACER in LDB1 
(open symbols) and LDB2 (solid symbols) (see lAve et all ( 20081 ) and lObermeier et al. 



(120111 )). Dashed lines indicate a power-law fit above 20 GeV amu 1 . 
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Fig. 



Boron-to-carbon abu ndance ratio as a fu nction of kinetic energy per nu 



cleon as measured 



3V TR ACER flObermeier et al.l l201lh . HEA P (lEngelmann et al.l ll990f). 
CR N dSwordv et al.lll99oh . ATIC JPanov et al.ll2007bh . CREAM fkhn et alJl2008h and AMS- 
01 (Uguilar et al.l boioh . Error bars are statistical (thin) and systematic (thick, only for 



TRACER). A simple model of the escape pathlength is indicated (dotted, see Eq. [7]). For 
the TRACER measurement the level of the subtracted contribution of atmospheric produc- 
tion of boron (dashed) is shown. 
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propagation index 5 



Fig. 3. — x 2 map in the parameter space of 5 vs. Ao for the Leaky-Box model fit to TRACER 
data. The best fit values are marked at (5, A ) = (0.53 ± 0.06, 0.311q3^ g cm -2 ) and the la 
contour is indicated. 
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Fig. 4. — Escape path length as a function of energy resulting from a fit to the boron-to- 
carbon data of TRACER. The dotted lines indicate the uncertainty range noted in Fig. [3J 
The dashed lines indicate the spallation path lengths of carbon and iron in the interstellar 
medium. 
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Fig. 5. — Abundance ratio of boron and carbon from TRACER (jobermeier et allboilh and 
from other measurements (see Fig. |2J). The solid line represents the best fit to the TRACER 
data alone; the dashed line is the best fit to all data combined. 



- 19 - 



E 

(0 



I — i i 1 1 1 1 



1 — i — i i 1 1 1 1 



1 — i — i i 1 1 1 1 



1 — i — i i 1 1 1 1 





o 

CM_ 



Spectrum at Earth 

Spectrum at Source 

■ TRACER 2006 
□ TRACER 2003 



^10 2 



in 

CD 



E 
> 

o 



in 
to 

ci 

LU 

X 

X 



10 



j, 



j i 



j i 



j i 



10' 



10 3 10 4 10- 

kinetic Energy [GeV/amu] 



Fig. 6. — Energy spectrum of oxygen as measured by TRACER. A source spectrum is fit to 
the data according to Eq. flHJ) with spectral index a = 2.37±0.12 (solid). The corresponding, 
observed spectrum is shown as a dashed line and uncertainty bounds are indicated as dotted 
lines. The lower uncertainty bound corresponds to a = 2.49, 5 = 0.59, and Ao = 0. The 
upper bound corresponds to a = 2.25, 5 = 0.47, and A = 0.8 g cm~ 2 . The normalization of 
the source spectrum is arbitrary. 



- 20 - 




i<r 1 1 10 10 2 io 3 10 4 

kinetic Energy [GeV/amu] 



Fig. 7. — Illustration of diffusive reacceleration models realized with GALPROP for different 
values for 5. The best model fit to all available cosmic-ray data has 5 = 0.34 (solid line). 



